Potential-driven eddy current in rippled graphene nanoribbons 
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It is well known that an eddy current will be induced in a conductor subject to a varying magnetic 
field. Here we propose another mechanism of generating nano-scale eddy current in rippled graphene 
nanoribbons(GNRs), which is only driven by an electric potential. In particular, it is found that 
under appropriate gate voltages, a local deformation may induce some unexpected global eddy 
currents, which form vortices in both rippled and entire flat areas of the GNR. We will explain that 
these vortices in flat areas is a manifestation of the nonlocality of quantum interference. 



A. Introduction 

Graphene, as the first experimentally available strictly 
two-dimensional crystal [1^ J , has offered a lot of new 
possibilities for both fundamental research and new tech- 
nologies due to its exceptional mechanical and electronic 
properties [3]-[5]. These single- atom-thick carbon mem- 
branes are very soft and flexible. Experiments have ob- 
served that graphene sheets are rippling even in free- 
standing conditionjBHS]. This feature is expected to have 
significant infiuences on their electronic properties and 
to result in some new observable effects[9-12 . Recently, 
Bao et al. reported a method to control and create nearly 
periodic ripples of sinusoidal form in graphene sheets [8^, 
which provides a possibility to design and manipulate 
electronic states in graphene sheets by controlling their 
local structure. According to the continuous model, a 
linear approximation of the tight-binding Hamiltonian 
in vicinity of its Dirac points, the electronic states of 
graphene can be described by a massless Dirac equation. 
In this model the action of a ripple can be represented 
by an effective vector field [T3UT5]. which essentially is 
the shift of the Dirac points due to the deformation, and 
a velocity tensor representing anisotropy of the energy 
band[16 . This continuous model is adequate if the de- 
formation is very smooth. However, this condition can- 
not be always satisfied in realistic graphene sheets, since 
some corrugations in graphene sheets are not so smooth. 

In this work, the local density of states(LDOS) and 
current fiow in rippled GNRs are investigated by us- 
ing the tight-binding model and non-equilibrium Green's 
function(NEFG) method[17-19 . It is found that the cur- 
rent distributions near each step edge of the conductance 
staircases are unstable, they will become non-potential or 
eddy flows if the GNRs are rippled, although their driven 
force is potential. In particular, there also occur some 
global eddy currents, which are still vortical in the fiat 
regions of the GNRs. The vortices in fiat areas is rather 
exotic from the classical view point, because there is no 
local defiection mechanism in these areas. We will ex- 
plain that this phenomenon is originated from a non-local 



quantum interference effect. This eddy current will give 
rise to a very inhomogeneous magnetic field in nanome- 
ter scale. It not only has crucial infiuences on the perfor- 
mance of graphene-based electronic devices and circuits, 
but may also have important potential applications. 



B. Model and Method 

The electronic properties of graphene can be described 
by the nearest tight-binding Hamiltonian 



H 



t(r-rOc'^(r)c(rO +H.C. 
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where r, denote two nearest lattice points, are elec- 
tron's creation and annihilation operators, t(r — r') is a 
distance-dependent hopping amplitude modeling the in- 
fiuence of deformation and can be fitted by 



t{v - v') ^ toe 



a(l-|r-r'|/0 



(2) 



where to = -2.75eV, a 3.37, / - 1.42A[20l|2T]. For the 
tight-binding model, we define the bond current j(r, r') 
from r to by the conservation equation r') = 

_d{n^t))_^ where n{rt) = 2c\rt)c{r't) is the electronic 
number operator, the factor 2 comes form the spin degree 
of freedom. From the equation of motion we have 



j{r,r' ,t) = -Re[t{r - r')G<{r%rt)], 



(3) 



where G^{Y't' ^rt) = i{c^ (Yt)c{Y't')) is the lesser Green's 
function. In steady states, the current can be written as 
j(r, r') = J jir^r' ^E)dE^ where 



J (r, v',E) = -Re[i(r - v')G< (r', r, E)] 



(4) 



is the current per unit energy [27. 

For the sake of simplicity, we only consider zigzag and 
armchair GNRs with a finite rippled region. The ripple 
is a static height fiuctuation given by 
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with (j)o an arbitrary constant and m an integer, the 
atomic in-plane displacement is ignored. A typical con- 
figuration of the rippled GNRs considered in this paper 
is shown in Figjl] 




FIG. 1: An N — 20 zigzag GNR with a sinusoidal ripple 
given by Eq.([5]), here /i :^ 2a with a 2. 46 A the length of the 
basis vectors, the ripple wavelength A = 30a, the included 
angle between the GNR axis and k is 30°, (j)o = — 7r/2, 
m = 1. 

In order to calculate the LDOS and current of GNRs 
by using the NEFG method [TTHTO]. the GNRs have to 
be considered as consisting of three connected parts: a 
rippled region as a sample and two semi-infinite ideal re- 
gions as leads [23, 24 . The left and right leads can be 
assumed to be in equilibrium states with different chem- 
ical potentials jj^l and jur^ respectively, where jul — J^r is 
very small. Thus the current can be written as 

j(r,rO = ^ rRe[t{r-r')G<{r\r,E)]dE, (6) 

where 

G< = G^S<G^ (7a) 

Qr,a ^ ^^±j _H^- (7b) 

with = E ± zO+ and He is the Hamiltonian matrix of 
the isolated sample. H^'^ = T^^£^ + S^'^ is the self-energy 
arisen from the couplings between the sample and two 
leads. ji = fL,R{^R L ~ l)' where the Fermi func- 
tions fL,R satisfy /lM^) - ^{i^l,r - E) at low tem- 
perature. S^^^ can be obtained by using the iterative 
method developed by Lopez-Sancho et.al.^. To reduce 
the calculation and memory requirement, we orderly de- 
note the layers of the sample from the left to the right by 
1,2, ••• thus the necessary sub-matrices G^^_i and 
^ for calculating current can be obtained by following 
back-and-forth recurrence procedure: 

=(^±/-ff„.„-S]^;j-i, (8a) 

^/?,n-l — ^n-l,n^n,n^n,n-l7 (8b) 

with n = l^ - ■ ■ ,2 and ^ = E^. For fiR < E < fiL, 

^n,n-l — 9n,nE^ri,n-lG^_l^^_l^ (9a) 



with n = 2, • • • , / and Gf ^ = iS<G? i, where 

GZ = [E^I - - S^i - Slfi)]-^ (10) 
Similarly, the LDOS 

p(r,^) = ^A(r,r), (11) 

where A = i{G'^ — G^)^ and the conductance 
2p2 

T{E) = -^Tr[r^(A - G^F^G^)], (12) 

with Tr^l = '^{^R L ~ ^R l)' ^^^^ be obtained by the 
same method. 

In order to see the ripple's influence on the current 
as clear as possible, the inherent deflection due to the 
lattice background has to be averaged out. To this end 
we deflne a cell-average current as an average vector of 
six bond currents around a honeycomb cell 

j(r.) = ig,(r.r,)^, (13) 

where Vij denote the vertices of the cell, Vc is the cell 
center. 



C. Results and Discussion 

First we compare the LDOSs and conductances for a 
flxed GNR with different ripples. Fig j2] shows the LDOSs 
and conductances for a zigzag and an armchair GNR in 
the presence of different ripples given by Eq.([5| with dif- 
ferent /I's and orientations (represented by slope angle 
6>k). The most remarkable change is that the conduc- 
tance near every step edge of the conductance staircase 
is remarkably decreased. Meanwhile, each corresponding 
van Hove peak in LDOSs is also broaden and succes- 
sively split into two sub-peaks when h reaches critical 
values, except for zigzag GNRs with 6>k = 0° (which will 
be further discussed in the last paragraph). The critical 
degrees of deformation can be represented by the associ- 
ated maximum relative bond elongation due to the ripple. 
This bond elongation can be roughly estimated from the 
ratio between the ripple's height h and wavelength A if we 
consider only the atomic height fluctuation. For zigzag 
GNRs with = 30°, this critical ratio is about 1/20, 
the corresponding maximum relative bond elongation is 
about 5%. 

To get an intuitive understanding of the microscopic 
behavior of electrons that gives rise to these changes, 
now we analyze the spatial distributions of the LDOS 
p(v^E) and current j(rc,^) at the corresponding ener- 
gies. Fig js] shows two examples of this kind of LDOS and 
current distributions in a rippled N = 20 zigzag GNR 
at O.SGeF and O.GleF, respectively corresponding to two 
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FIG. 2: (a) LDOS and (b) conductance of rippled AT = 20 
zigzag GNRs as shown in Fig[l] for different /I's, where A = 
30a, = 30°. (c) LDOS and (d) conductance of identical 
GNRs as Fig[l]for different ^k's, where h — 1.5a, A = 30a. 
(e) LDOS and (f) conductance of A/" = 20 armchair GNRs 
for different /I's, here A = 40a, 6'k = 20°. (g) LDOS and 
(h) conductance of identical GNRs as (e,f) for different ^k's, 
where h — 1.5a, 00 — — 7r/2. 

sub-peaks split off from the first van Hove peak in con- 
ducting sub-band. We can see from the LDOS distribu- 
tions (Figj3ja,c)) that the lower sub-peak is mainly local- 
ized in the ripple area; while the upper one (belonging to 
the bottom of the first conducting sub-band) is extended, 
its LDOS does not decay outside the ripple region. More 
importantly, their current distributions both occur re- 
markable vortices (Figj3jb,d)). Some current lines even 
form closed loops, showing that these currents are not 
potential flow. Similarly, the second van Hove peak also 
splits into two sub-peaks and their current distributions 
also have this vortical feature. The valence bands also 
occur identical phenomena. When the bias is reversed, 
these eddy currents will also be exactly reversed, thus an 
alternative bias can drive a varying eddy current in this 
kind of rippled GNRs. 

At first sight, these vortices seems can be ascribed 
to Landau-like quasi-bound states caused by the ripple- 
induced pseudo-magnetic field, since according to the 
continuous model a geometrical deformation will induce 
an effective vector potential 

^(V3(t3-t2),t3+t2-2ti) (14) 

with ti^2,3 three nearest hopping amplitudes [15 . Actu- 
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FIG. 3: (a) LDOS and (b) current distribution at 0.56ey (the 
lower sub-peak split off from 1st van Hove peak) in an A/" = 20 
zigzag GNR as shown in Fig[l] (c) LDOS and (d) current at 
0.6 ley (the upper sub-peak of the 1st van Hove peak) in the 
same GNR. The brightness in (a,c) represents the LDOS, the 
contours depict the ripple. The brightness in (b,d) represents 
the height of the ripple. 



ally, some conductance fluctuation phenomena similar 
to the results shown in Fig.Q have been observed in 
the measurement of differential conductance of graphene 
sheet with nanobubbles, and the researchers interpreted 
all this conductance fluctuation as the contribution of the 
Landau levels arising from the pseudo-magnetic field [TJ. 
It is true that most of the above results can be explained 
by this pseudo-magnetic field. For example, if k is exactly 
along some special directions, e.g., = 0°,30°,60°, 
we can deduce some general properties of A by quali- 
tative analysis, which can be compared with the numer- 
ical results. In these special cases the A will has the 
form Aq cos(k • r + (/)o) and will be parallel or perpen- 
dicular to k[16 . More specifically, for rippled zigzag 
GNRs of l9k = 0°,60° or armchair GNRs of l9k = 30°, 
A is nearly parallel to k, so the pseudo-magnetic field 
V X A c:^ — k X Aq sin(k • r + (/)o) is very small, thus 
the ripple influence will be much weaker than other rip- 
ple orientations; while for armchair GNRs of ^k = 0° or 
zigzag GNRs of ^k = 30°, A is perpendicular to k, so the 
pseudo-magnetic field is strong and will strongly disturb 
the electronic states. For these special cases, the above 



4 



explanation is qualitatively in accord with the result in 
FigJIIthere are still some more subtle problems, which 
will be discussed in last paragraph). 

However, these current distributions also exhibit a re- 
markable feature revealing that this is not the whole 
story. We notice that if the energy slightly higher than 
the bottom of a conducting sub-band (or lower than the 
top of a valance sub-band), the vortices of the eddy cur- 
rent will not entirely lie within the ripple region, but 
also occur in flat regions far from the ripple, where the 
pseudo-magnetic field has vanished, as can be seen from 
Figjsjd). This feature reveals that these vortices can- 
not be interpreted as Landau states produced by the 
pseudo-magnetic field, because all ripple-induced scatter- 
ing mechanisms, both the pseudo-magnetic field and the 
velocity variation, only act locally within the ripple area, 
hence the electronic waves should freely propagate in fiat 
areas. Figj4ja,b) give other two examples of this global 
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that both of them have pronounced vortices in the fiat 
areas far beyond the ripple, where the pseudo-magnetic 
field has certainly vanished. Generally, if the ripple is 
slope relative to the GNR axis, the current distributions 
within the energy range slightly above the bottom of a 
conducting sub-band or below the top of a valence sub- 
band will occur remarkable vortices in entire fiat areas. 
These vortices appeared in fiat areas are rather exotic be- 
cause there is no local responsible defiection mechanism. 

In order to explain the origin of these exotic vortices, 
we have to notice two basic properties of the electronic 
states of GNRs under perturbation. The first is that 
the electronic states will become superposition of par- 
tial waves with approximate energies. The second is that 
the velocity direction (forward or backward relative to 
GNR's axis) of the states near the bottom of a conduct- 
ing sub-band (or the top of a valence sub-band) is unsta- 
ble, because its velocity dE{'k)/dkx {kx is the momentum 
component along the GNR axis) is very small, so its sign 
can be easily changed by a small variation of E(k), as 
illustrated in Figj5|a,b). Therefore, if there is a nonuni- 
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FIG. 4: (a) Eddy current in the same GNR as Fig[T]at leV, 
the upper sub-peak of the second van Hove peak, and (b) 
1.37eV, the upper sub-peak of the third van Hove peak, (c) 
Eddy current at O.GeV in a similar GNR as Fi g|3[c,d) but 
with a hyperbolic-surface ripple given by Eq.(15). (d) LDOS 



and (e) current distribution at 0.69eV in an N — 20 armchair 
GNR with a ripple given by Eq.([5]), where /i = 1.5a, A = 40a, 
= 30°, 00 = — 7r/2. (f) LDOS and (g) current distribution 
at 0.631ey in a GNR as (d,e) except 6'k = 0°. 

eddying currents in the same zigzag GNR as FigjS) their 
energies are slightly above the second and third van Hove 
peaks of the conducting band, respectively. We can see 



FIG. 5: (a) Conducting band of graphene near a Fermi point, 
the red dashed line represents a sub-band of a GNR. The ve- 
locity of states near the minimum of the sub-band (red dashed 
line) can be easily reversed under a small variation of the en- 
ergy band, as shown in (b). (b) Contour map of energy bands 
for a perfect (blue-purple) and a deformed (yellow) graphene. 
For the deformed one, only two bonds are elongated respec- 
tively by 5.3% and 2.8%. The red(yellow) arrow represents 
the velocities of a state in perfect (deformed) graphene. (c) 
Current distribution of three plane waves. 

form deformation, the perturbed state corresponding to 
an originally forward state in these energy ranges may 
contribute a backward fiow in the deformed regions. In 
particular, this local backward fiow demands the per- 
turbed state to include backward partial waves in fiat 
areas in order to satisfy continuity condition at the in- 
terfaces between deformed and fiat areas, though there 
exists no local responsible scattering mechanism. These 
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partial waves of different directions are always spatially 
superimposed due to the edge reflection and will produce 
interference. Consequently, the current density no longer 
equals to the sum of the currents of every partial waves, 
j(r) 7^ ^(k|j(r)|k), but must include interference terms 
(k|j(r)|k'), similar to the usual interference for proba- 
bility density. These kind of current interference pat- 
terns will be winding or eddying flows in areas without 
any local deflection mechanism. As a simplest example, 
Figjsjc) shows the current interference pattern of three 
plane waves ?/^(r) = e^^^'^ + g*k2 r _^ c^^s-r^ where ki^2,3 
have equal length while their included angles are 120°. 
We can see that its current distribution forms a very sim- 
ilar eddying pattern. Obviously, the vortex scale of these 
eddy currents arising from the interference is in propor- 
tional to the wavelengths of the partial waves. For the 
eddy currents in rippled GNRs, this character can also be 
verified by comparing Figj3];a,d) and Fig|4];a,b). We find 
that the vortices become smaller and smaller with the in- 
creasing of energy owing to the linear dispersion relation 
E (X k. Conversely, the vortex scale will be very large in 
wider GNRs because the energies of each corresponding 
step of the conductance staircases will be smaller, it is 
only limited by the electronic interference length. 

According to this explanation, these global eddy cur- 
rents would be a rather ubiquitous effect in rippled 
GNRs, although their patters depend on specific rip- 
ple configurations. Figjljc) is another eddy current in 
an identical zigzag GNR in the same energy range as 
Figjsjb) but with different ripples, which is a hyperbolic 
surface 



zir) 



h[l-f{v)], /(r)<l 
0, /(r) > 1, 



(15) 



where/(r) = ^ (^) + (^) , with h = 4a, ao = 

15a, 60 = 25a, §1,2 = (V2/2, =FV^/2)(its maximum bond 
elongation is about 3%). By comparing Fig jsjd) and 
Fig|4jc), we can see that their current distributions in 
flat regions are very similar although their ripple are 
very different. Similar to this structural insensitivity of 
the global eddying character, it is conceivable from the 
above explanation that this vortical character will also 
not very sensitive to the energy variation. Actually, the 
representation of these eddying states in the LDOS curve 
is not a sharp peak like quasi-bound states, but a broad 
and smooth one forming a piece of continuous spectrum 
(Fig§. 



In addition, there are few special cases worth to be 
particularly pointed out. The first case is zigzag GNRs 
with = 0°- In this case the current lines remain to 
be straight lines and no vortex occurs, because in this 
case the A is parallel to k according Eq.(14) or [16], so 



the V X A can be ignored; moreover, the incident cur- 
rent is along a symmetric axis of the energy band, so 
the refraction also does not change its direction. The 
second is zigzag GNRs with Oi^ = ±60° or armchairs 
GNRs with = ±30°. Similar to the first case, here the 
V X A are also very small, however, there exist appar- 
ent vortices, as shown in Figjlj^e), because in these cases 
the anisotropic deformation of the energy contours (see 
Figjsjb)) will result in similar backward flow. The third 
one is armchair GNRs of = 0°, although the ripple 
induce a strong pseudo-magnetic field and results in ap- 
parent quasi-bound states (Figj2|f) and FigQf)), but 
their currents do not occur any vortex(Fig|4]^gjy, because 
the direction of the incident current is the symmetric 
axis of two pseudo-magnetic fields for two inequivalent 
Dirac points (mutually symmetric) as well as the prin- 
ciple axis of the velocity tensor, so the action of the 
pseudo-magnetic fields of two Dirac points will be mutu- 
ally canceled out and the refraction and dispersion due to 
the velocity anisotropy also does not change the current 
direction. 

In summary, the current flows near every step edge of 
the conductance staircases of a GNR are rather unsta- 
ble. They will become eddy currents if there occurs a 
slope ripple. These eddy currents can be divided into 
two classes. The first one are carried by Landau-like 
states caused by the pseudo-magnetic field, these states 
have slightly lower energies and their current distribu- 
tions form vortices only within the ripple region. In con- 
trast, the second one are carried by some special scatter- 
ing states, which have slightly higher energies and include 
backward partial waves in flat areas. Consequently, they 
will form some global eddy currents due to the interfer- 
ence of these partial waves. This global eddy current is a 
manifestation of the non-locality of quantum interference 
effect. 
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